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Self-propelled motion of a fluid droplet under chemical reaction 

S. Yabunaka.^y T. Ohta.^y and N. Yoshinaga^f)] 

^''Department of Physics, Kyoto University, Kyoto 606-8502, 

Japan 

^'^WPI-AIMR, Tohoku University, Sendai, 980-8511, Japan 

(Dated: 24 February 2013) 

We study self-propelled dynamics of a droplet due to a Marangoni effect and chemical 
reactions in a binary fluid with a dilute third component of chemical product which 
affects the interfacial energy of a droplet. The equation for the migration velocity of 
the center of mass of a droplet is derived in the limit of an infinitesimally thin inter- 
face. We found that there is a bifurcation from a motionless state to a propagating 



n , I state of droplet by changing the strength of the Marangoni effect. 
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I. INTRODUCTION 

Self-propelled motion of particles has attracted much attention recently from the view- 
point of non-linear physics far from equilibrium. There are several experiments of self- 
propulsion of droplets in fluids^"-. It has been shown that the Belousov-Zhabotinsky reac- 
tion composed in a fluid droplet triggers a spontaneous motion of a droplet^. Computer 
simulations of convective droplet motion^ and nano-dimer motors^"^ driven by chemical re- 
actions have also been carried out. There are theoretical studies of droplet motion due to 
an interfacial tension gradient along the droplet surface^i^i^i^. However, these theories are 
concerned only with the steady velocity of a droplet. As a related theoretical study, the 
mesoscopic description of the thermo-capillary effect has been formulatec^ii. A transition 
between a motionless and migrating droplet driven by chemical reactions has been studied 
in a system where a droplet is on a solid substrate^^. 

It should be noted that self-propelled motion of particles has been investigated in a 
different field of physics. It has been known that a pulse or a domain in excitable reaction 
diffusion systems exhibits a bifurcation from a motionless state to a propagation state by 
changing the system paTametersM-^. A reaction-diffusion system is represented by a set of 
nonlinear partial differential equations, that is often investigated by numerical simulations 
due to the limitation of analytical calculations. Nevertheless, the theory of domain dynamics 
in the vicinity of this drift bifurcation has been developed, e.g., for the interaction between 
domains^""— and for deformations of domain^"—. 

The purpose of the present paper is to extend the previous studies in reaction-diffusion 
systems to the droplet motion in chemically reacting fluids. We introduce a model system of 
binary fluids where a chemical reaction takes place inside a droplet. The chemical component 
produced diffuses away from the droplet and influences the interfacial energy. The long 
range hydrodynamic effects are treated with a Stokes approximation supposing that the 
relaxation of the fluid velocity field is much faster than that of the concentrations and that 
the Reynold number is sufficiently small in the system considered. We will show that there is 
a drift bifurcation at certain threshold of the Marangoni strength as in the reaction-diffusion 
systems mentioned above. The time-evolution equation of the center of mass of droplet is 
derived near the drift bifurcation by taking into consideration of the hydrodynamic effects. 

In the next section (section [Tl|) , we describe our model system and the interface dynamics. 



The equation of motion for the center of mass is derived in section IIIII Discussion is given 
in section IIVI The force acting on the droplet interface is formulated in Appendix \M Some 
of the details in the derivation of the velocity of the center of mass are given in Appendix 
[Bl The formulas used in the evaluation of the coefficients in the time-evolution equation 
for a droplet are summarized in Appendix O The convective effect of the third chemical 
component is estimated in Appendix 

II. MODEL AND INTERFACE DYNAMICS 

We consider a fluid mixture where the free energy is given in terms of the local concen- 
tration difference = 0^ — 0b by 



F{0} = / dr 



■(V0)^ + /gl(0) + /o(c) 



(1) 
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where (I)a{4'b) is the local concentration of the component A (B) and /o(c) = cine. The 
coefficient i? > is assumed to depend on c as B{c) = Bq + Bic with Bq and Bi constants 
and /gl(0) is a function of such that phase separation takes place at low temperatures. 
Here we have assumed existence of a dilute third component whose concentration is denoted 
by c. The logarithmic term (/o(c) = cine) arises from the translational entropy of the dilute 
component. The spatial variation of c is also assumed to be broad enough compared to that 
of which constitutes a sharp interface. 

The time-evolution equation for is given by 

_| + V.W) = V^-, (2) 

where v is the local velocity whose equation is given by eq. (jl]) below. Hereafter we consider 
an isolated droplet such that the concentration variation is 0(x) = 0e > inside the droplet 
and 0(x) = — 0e at the surrounding matrix. The equilibrium value 0e is determined by 
equating the rhs of eq. (E]) to zero. The dilute component c is assumed to obey 

_ + V ■ (t;c) = DV^c - 7(c - c^) + Ae{R- \r - rbl) , (3) 

where 9{x) is the step function such that 9{x) = for x < and 6{x) = 1 for x > 0. The 
first term on the rhs of eq. ([3]) arises from V ■ [L{c)'V6F/5c] with L(c) = Dc where D is 
positive constant. The c-dependence of the Onsager coefficient L is necessary for a dilute 
component^i. The second term in eq. ([3]) indicates consumption of c with the rate 7 > 



due to a chemical reaction and with c = Coo for \f\ — )■ oo whereas the last term represents 

production of c, which occurs inside a droplet with radius R, whose center of mass is located 

at fc- In the most parts of the present paper, the coefficient A is assumed to be positive 

and stands for the strength of the production. However the theory can also hold for A < 

with a slight modification. 

The Stokes approximation is employed for the local velocity v and it takes the form 

-^SF ^6F 
= - Vp - 0V-— - cV-- + m^^v, (4) 

d(p oc 

where p is determined such that the velocity field satisfies the incompressibility condition 

V ■ V = 0. The viscosity tjq is assumed, for simplicity, to be a constant independent of 0. 

The force arising from the first, second and third terms can be written as 

SF 5F 

p = _v> - 0V"— - cV"— = -W + /|f + fl, (5) 

where p" has some additive terms to p, whose explicit form is unnecessary for incompressible 
fiuids since only the transverse components of the velocity is relevant. In Appendix |Al we 
show that the normal and tangential forces are given, respectively, by 

/|f = -n"5(c)|V0|2(V-n), (6) 

/! = (5./3-n"n'3)(V^5)|V0p, (7) 

where the unit vector n is directed to the outside of the droplet, i.e., n = — V0/|V0|. 
The repeated indices imply the summation. When we are concerned with the large scale 
compared with the interface width (or the sharp interface limit), the factor |V0|^ is localized 
in the interface region. In this situation, the forces are localized on the interface at a which 
denotes a location on the interface so that we may rewrite Eqs. ([6]) and ([7]), respectively, as 

/|f(a) = -nV(a)(V-n), (8) 

/^(a) = (5„^-n"n/^)(VM/, (9) 

The interfacial tension is defined by 

aia) = JdwBic) ^^J ^ B,{a) j dw (^^\ (10) 

where w is the coordinate along the normal to the interface and Bj is the value of B at 
the interface. It should be noted that the derivative in V^a is not restricted to the two- 
dimensional space on the interface regarding cr(a) as cr(c(r)). After taking the derivative in 



three dimensions, we may take the value on the interface. This interpretation is consistent 
with Eq. ([7]) in which V'^ acts on the weak spatial variation of c. The tangential component is 
automatically extracted by the projection {dap — 'nf^n^)- Equations ([H]) and ^ are consistent 
with the boundary condition employed in hydrodynamics with multi-component fluids^^. 

Substituting Eq. (|5]) into Eq. (j4]) and using the incompressibility condition, the local 
velocity of fluid is given by 



y"(f,t) = / da'T""^ {r,r{a!))n^{a')a{a\t)K{a\t) 

da'T^^ {f,r{a')) [5p^ - n^{a')n\a')]{V^a)i, (11) 



where da' is the infinitesimal area on the interface. The integral is taken all over the interface. 
The Oseen tensor is given by 



T"'^ ff, r" 



1 r. s'^sP 



(12) 



QnrjQS 

with s = f — r* . The mean curvature is defined hy K = —V ■ n. 

The right hand side in the time-evolution equation ([2]) for can be ignored when the 
hydrodynamic effects are dominant^^. From the left hand side of Eq. ([2]), we note that the 
normal component V{a,t) of the interface velocity is given by 

\/(a,t) =t;"(f(a),t)n"(a). (13) 

Substituting Eq. ( ITTj) into Eq. ( iT3|) . we obtain 

V{a,t) = V^ + V2, (14) 



where 



Vi= da'n"{a)T"^ {f{a), f{a')) n^(a>(a', t)K{a', t), (15) 



and 



Vo 



f da'n''{a)T''^ {r{a),ria')) [5p^ - n^(a>^(a')](V^(T)/. (16) 



The velocity u of the center of mass of an isolated droplet can be obtained from V{a, t). 
The geometrical consideration leads te^ 



u'' = ^ I daV{a)R'^{a), (17) 



— * 

where fl is the volume of the droplet and R{a) is the position vector directed from the center 
of mass to the interface. For a spherical droplet with radius R, we have Q = AnR^/S and 
R{a) = n{a)R. 

In order to determine the migration velocity u, we have to evaluate the interfacial tension 
and its spatial derivative as Eqs. (fT5|) and ( !T6|) . which may depend on the concentration c. 
In this way, we take into account the Marangoni effect. To this end, we assume that the 
interfacial tension depends on c/ as 

a = (To + fxic/, (18) 

where ctq and cxi are constants determined from the expression oi B = Bq + Bic. However, 
the explicit form of (Tq and ai as a function of Bq and Bi are unnecessary in the argument 
below. Substituting (TT8l) into ( TT7|) . we obtain for a spherical droplet with K = —2/R 

u" = u'^ + u^, (19) 



where 



u 



2cti 



dan'^ia) da'n\a)T^\r{a),r{a'))n\a')ci{a'), (20) 



u° = '^ I dada'n''{a)n\a)T^^ {r{a),r{a')) 



n 

X {6p^-n^{a')n^a'))\/^c{a'). (21) 



Equations (!20|) and (12T|) are derived in Appendix [B] as 



< = -^Y^ / da'n"ia')cjia'), (22) 



.-^/<'"'(*».-"°M.V))Kc),. (23) 

In the next section, we will derive the time-evolution equation for u from Eq. (lT9l) with ( l20ll 
and ( pT]) by solving Eq. (|3]) for the third component c. 

It is remarked that, when c{r) is set as c = cq + Ciz instead of solving Eq. (|3]), we obtain 
from Eq. (IT9l) with (l20l) and (12T1) the stationary migration velocity n = — 2(TiCii?/(15ryo) 
which agrees with the known result obtained by the conventional theory of the Marangoni 
effect^. 
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III. EQUATION OF MOTION FOR A DROPLET 

In this section, we derive the equation of motion for a droplet. Since the major hydro- 
dynamic effects have been taken into account as in Eqs. flT^ . flT^ and (1161] . we ignore the 
convective term V ■ (vc) in Eq. 1^. We will show in Appendix [D] and in section HVl that this 
term causes a shift of the bifurcation threshold but is not expected to change the bifurcation 
behavior essentially. 

The configuration of the component c around a droplet can be obtained by solving the 

following equation 

oc 

— = DV'c - 7(c - Coo) + Ae{R- |f - rbl) . (24) 

Hereafter, we consider the case oi A > that the component c is produced inside the droplet, 
diffuses away, and vanishes at |r| — )■ oo i.e., Coo = 0. The method can also be applied for 
A < with the boundary condition c = Coo 7^ for |r| — )■ 00. In terms of the Fourier 
transform, Eq. (12^ can be written as 

^ = -D{q' + ^')c,+ H,, (25) 

where 

^=S)*- P*" 

H,= AS,e'''<--, (27) 

with the form factor of a sphere 

S^^= f d^rexp{tq-r)9{\f\-R) (28) 

^ ^^ sm{qR) - qRcos{qR) 

The Fourier component c^ has been defined as 

Cg-.= /dVc(f)e^«^^ (30) 

By assuming the relaxation of the composition c is sufficiently rapid compared to the 
motion of interface, we solve Eq. (l25|) by means of an expansion in terms of the time 
derivative. 



JO) , ^(1) , (2) , fS) 



'^'+c^' + c^'+^' + , (31) 



Q q q q 



where we have defined 

The short time expansion (13T!) is justified in the vicinity of the supercritical drift bifurcation 
where the velocity of a droplet n = |n| is arbitrarily small. That is, the smallness parameter 
of this expansion is given by 

e = ^«l. (33) 

where the denominator is the characteristic time of c. After the inverse Fourier transform, 
the composition c/ at the interface is given by 

ci = cf [re + ^ + cf (r-G + s) + cf (fc + s) + cf (rb + s) , (34) 

where 

cf{fG + s) = A fGgSge''^-''^e-''^<''^^^ = A fCgSge-''^', (35) 

J q J q 

c^i\rG + s) = -A jjzq- n)G^5,e-^^^-^^ = u'^^Q^is), (36) 

cf\rG + ^ = A fpq- u)GlSge-''^-' + A f{tq- ufGlSge-"^-' 

J q J q 

= -u-^Qs{s) + u-u^^^^Q,is), (37) 

cf\fG + s) = -A f{zq- ufG^^Sge-^'^-' 

Jq 

The terms with the higher order time derivatives have been ignored. The migration velocity 
is given by 



u = — . (39) 



We have defined Qn{s) by 



Qn{s) = A / G^^Sge-'^'^^ (40) 

Jq 



Since we have obtained the concentration profile of c for a given interface configuration, 
we can now evaluate the velocities in Eqs. (!22|) and ( 123|) . which are carried out in Appendix 
[Cl It turns out that there is a simple relation U2 = — (3/4)ui. From the results obtained 
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FIG. 1. The scaled coefficient m as a function of R. 
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FIG. 2. The scaled coefficient f as a function of R. 

in Appendix [Cl the time-evolution equation for the center of mass is given up to the cubic 
non-hnearity by 

|2 



where 



m = -M 
T= -M 



mil"' = (—1 + r) u" — gu"' \u\ 
dQs 



ds 
dQ2 



ds 
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3M 



s=R 

1 
s=R 

2 dQi 2 d'^Q^ d^Qi 



with 



i?^ ds R ds"^ ds^ 
2ai 



s=R 



M = 



15r]o' 



(41) 

(42) 
(43) 
(44) 

(45) 



As will be shown below, all the coefficients m, r and g are positive. The term proportional 
to u'^ does not appear, because it is not a dissipative term. The third order term — (? |u| u" 



0.02 
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FIG. 3. The scaled coefficient ^ as a function of R. 

is needed to make the migration velocity finite. By choosing 1//3 as the characteristic length 
and 1/{DP^) as the characteristic time of the problem, Eq. (14T]) can be written in terms of 
the dimensionless quantities as 



m- 



di 



-Tc + f)u 



gu \u\ , 



(46) 



where i = tD(3'^, u = u/{D[3) and 



Tc 



D'/3^ 15rioD'/3 



2/33 



MA 



2aiA 



(47) 



Here we consider the case that aiA is positive. It is remarkable that all the parameters in 
the system are combined together as r^ given by (1471) so that Tc is the only dimensionless 
parameter. This is the case even if one takes account of the convective term in Eq. (I3l) since 
it does not contain any extra parameters. The dimensionless coefficients depend only on 
R = up and are given by 



m{R) = mD(3^Tc, 

f{R) = TTc, 

g{R) = g{Df3)\. 



(48) 
(49) 
(50) 



These scaled coefficients have been evaluated numerically and plotted in Figs. [H [2] and [3] 
which indicate that those are definitely positive. 
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FIG. 4. Translational motion of a droplet. The droplet is migrating to the right under the non- 
uniform distribution of the c component indicated by the small dots. 

IV. DISCUSSION 

We have formulated the theory of self-propulsion of a droplet caused by a Marangoni 
effect and chemical reactions. Equation of motion for a spherical droplet has been derived 
as Eq. ( I46l) which exhibits a drift bifurcation. The hydrodynamic effects are taken into 
consideration by the Stokes approximation for the fluid velocity. This is justified when the 
time variation of the concentrations is much slower than that of the local fluid velocity. We 
have made two assumptions. One is the assumption that the interface (surface of droplet) 
is infinitesimally thin. This assumption is satisfied when the droplet radius is much larger 
than the interface width. The other assumption is that the relaxation of the component c is 
much faster than the interface motion. Since the interface velocity is arbitrarily small in the 
vicinity of the drift bifurcation threshold, the second assumption is consistently justified in 
the theory. 

The mechanism that a droplet undergoes a translational motion in our model for A > 
and 0"! > is as follows. When a droplet is motionless, there is an isotropic concentration 
distribution of c around it. The concentration profile outside the droplet is a decreasing 
function of the distance from the center of mass. Let us suppose that the position of the 
droplet is shifted slightly. Then, the concentration of c decreases (increases) at the front 
(rear). If the relaxation rate of the component c is infinite, this concentration unbalance is 
recovered instantaneously. However, when the relaxation is finite, the droplet tends to shift 
further since the interfacial energy is an increasing function of c. This is shown schematically 
in Fig. m In fact, it is found that the terms with the coefficients r, m and g in Eq. ( 1411) arise 
from the higher order terms (c^ , c^ , and c^ , respectively) in the short time expansion 
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in Eq. (13T]) . Therefore, if the time-delayed effect ru" dominates the term —u" which 
corresponds to the Stokes drag force, the droplet undergoes migration. It is noted that this 
argument can also be applied to the case A < and cti < 0. 

We can estimate the effect of the convective term in Eq. ([3]) which have been ignored in 
the treatment in section IIIII In Appendix [Dl we derive the correction from the convective 
term up to the first order of the perturbation expansion. The coefficient r is evaluated since 
this quantity is directly related to the drift instability threshold. In the limit i? — > 0, we 
obtain 

r = i^P. (51) 

Tc 15 ^ ' 

When the convective term is not considered, we have P = 1 from (jSj). The first order 
correction from the convective term gives us P = 31/56 as shown in Appendix [Dl Since 
migration of droplet occurs for r > 1, this indicates that the stronger Marangoni effect is 
necessary when the convection of the third component exists. 

The reason as to why the convective term of H{r) = v-Vc tends to suppress the Marangoni 
effect can be understood as follow. Substituting the local velocity given by Eq. ( ]D2[) . we 
have the value of H at the interface 



Hi = u-S/c 



I 



(52) 



When A is positive, Vc/ and u are anti-parallel (parallel) to each other at the front (rear) of 
the moving droplet so that we may expect that if < (>)0 at the front (rear) area. Since the 
first order correction to the concentration c is given by c{r) = —[—DV^ + ^]~^H{f) and the 
operator [— Z^V^ +7]~^ is positive definite, the concentration c tends to increase (decrease) 
at the front (rear). This is just opposite to the concentration variation described above for 
the mechanism of translational motion. 

One of the characteristic features of the present theory is that all the parameters in the 
model equations are combined as r^ given by Eq. (HTj) which determines the threshold of the 
drift bifurcation. Since Tc is inversely proportional to A and ui, the self-propulsion is easier 
for the stronger production of c (i.e., larger values of A) and for stronger Marangoni effect 
(i.e., larger values of ui). Note that f is an increasing function of the radius of droplet. This 
means that the drift instability is favorable for larger droplet if other parameters are fixed 
and if any shape instability would not occur. 
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We make a remark on the sign of the Marangoni factor. We have restricted ourselves to 
the case of Aai > 0. When this quantity is negative, the coefficients m and g are negative in 
Eq. fl4ip . Therefore, in this case, we have to take account of the higher time derivatives and 
the higher nonhnear terms of u. However, this is beyond our present theoretical formulation. 

In the present theory, the third component is produced inside a droplet. However, if it is 
produced only on the droplet surface, the step function in Eq. ([3]) should be replaced by the 
delta function. We expect that the results obtained in the present paper are not essentially 
altered if the component c diffuses to the inside of droplet as well as the outside. Such a 
model has been studied where the time-evolution equation of surfactant on the surface of 
droplet is introduced explicitly^^^. 

A self-propulsion of an oily droplet has been observed in a micron size^. In this ex- 
periment, the molecules which constitute the droplet are produced by a chemical reaction 
which takes place at the droplet surface. Another experiment by Thutupalli et alA shows 
that an aqueous droplet of the order of lOO/xm surrounded by oil with surfactant molecules 
undergoes migration by causing a non-uniform surface tension due to bromination on its 
surface. In these experiments, however, it seems that the bifurcation from a stationary state 
to a moving state predicted in the present study has not been observed. Further systematic 
experiments are desired. 

Since fluid droplets are soft, they are generally deformed in migration. A coupling between 
migration velocity and shape deformations has been formulated recently in an excitable 
react ion- diffusion systent^. Extension of such a theory to the present hydrodynamical sys- 
tem will be carried out in the future. 
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Appendix A: Derivation of the forces 

In this Appendix, we derive the formulas (|6]) and ([7]). The force ([5]) is written as 

6F 6F 

/° = -V>-0V°---cV"--. (Al) 

0(f) oc 

Substituting the free energy ([1]) into Eq. flAll) . we obtain the modified pressure 

P =P + (f^-Q— - /gl + c-^ - /o, (A2) 

and 

/" = -VV - (V°0)(V»(V^5) - 5(c)(V°0)VV + ^(V"5)|V0|' 
= -W + \V<P\\V^B){6^p - n'^n^) 
- Bic){V''<P)V'<P - ^(V"i?)|V0r, (A3) 

where n = — V(/)/|V0|. In the last term on the first line of Eq. (1A3I) . we have used the 
relation {V"'c)dB/dc = V"-B. Note the formula 

VV = -V^ (ri^|V0|) = -(V • n)\V(j)\ + nf^V^\V(j)\ 
= -(V ■ n)\V(j)\ + n''V^n^V^(j)) 

= -(V-n)|V0|+nV(V^V'^0), (A4) 

where we have used the fact that n^CV^n^) = (l/2)V'^(n'^)^ = since (n^)^ = 1. Substitut- 
ing this into Eq. ( 1A3I) . we obtain 

fa = _VV + |V0p(V^5)((5,^ - n"n^) - 5(c)n"(V ■ n)|V0p 

- -(V"5)(V0)2 + fi(c)n"|V0|n%'^(VW^0) 

+ Iv" (5(V0)2) - iv" (5(V0)^) 

= -VV' + |V(/.p(V^5)((5„;3 - n"n^) - 5(c)n"(V ■ n)\V(p\^ 

- B{c)n''\V(j)\n^n^{WV'^<j)) + B{c){V''V^(l)){V^(j)) 

= - W + |V</.p(V^5)(5,;3 - n'^n^) - B{c)n''{V ■ r2)|V0|' 

+ S(c)(VW^0)(V^0)(5,^ - n°n^), (A5) 



where 



p" = p> + ^B{Vc^f. (A6) 
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Therefore the force / can be divided into the normal and the perpendicular components 

r = -VV + /|f + /±, (AT) 

where 

/|f = -n"S(c)|V0|2(V-n), (A8) 

/^ = (5a/3-n"n^)[(V^i?)|W|2-5(c)|V0|(V^V^</))n^]. (A9) 

The second term in Eq. (]A9P is negligible compared to the first term in the sharp interface 
limit. In fact, we have 



|V0|(V^V^0)n^] = |V0| {V^n^)\V(f)\+n^{V^\V(j)\) 
1 



7/3|V7A|2 



nP' 



-V'\V<P\\ (AlO) 

where we have again used the formula n'^(V'^n'') = (l/2)V^(n''')^ = 0. The integral of 
-B(c)V'^| V0p across the interface vanishes provided that B varies weakly across the interface. 
Therefore we ignore the second term in Eq. (jAOp . 

Appendix B: Derivation of the migration velocity 

In this Appendix, we derive Eqs. ( 12 2 p and (1231) . In order to obtain Eq. (I22p . the following 
formula for a spherical droplet^^ is necessary. 

j da'n''{a)T''^{r{a),r{a'))n\a')YUa') = EiYUa). (Bl) 

where 

' r/o(2/-l)(2/ + l)(2/ + 3)' ^ ' 

and Yjm(a') is the spherical harmonics. The representation of the unit vector n in terms of 
Yim is also necessary. 

n = (sin 6 cos </>, sin 6 sin </>, cos 6) (B3) 

-(-Fn + ri_0,yy(rn + Fi-i),yyFioj. (B4) 

Applying these formulas to Eq. (!20|) . one can carry out the integral over a so that Eq. (122|) 
is obtained. 
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Next we calculate Eq. (12T]) . First we make an ansatz as 



da'T""^ {f{a),r{a'))n^{a')n'{a') = X5^^ + Yn''{a)n^{a 



(B5) 



The unknown constants X and Y are determined as follows. We note the identities; 

da'T"^ {r{a),f{a')) n^(a')n°(a') = 3X + Y, (B6) 

dan"{a)T'"^ {f{a),r{a)) n''(a')n^(a)n^(a') = X + Y. (B7) 



/' 



The left hand side of these expressions is readily evaluated as 

da'T'"^ {f{a),f{a')) n^(a')n"(a') 
l + sm'^{9/2) 2R 



D /"I 

-— / d{cose)- . 

8t]o J„i sm(6'/2) 



3?7o' 



(B8) 



da'n''{a)T''^ {f{a),f{a')) n^{a')n\a)n^{a') 

SvoJ^i ^ ^ sin(^/2) 15r/o' 

where 9{> 0) is the angle between n[a) and n[a'). Therefore we obtain 



X 



Y 



R 

5rjo' 
R 

15?7o 



By using the formula (IBSp . Eq. ( 123|) is readily obtained. 



Appendix C: Derivation of the coefficients 



(B9) 

(BIO) 
(Bll) 



In this section, we derive the migration velocities by evaluating Eqs. 
Substituting Eqs. (jSS]), ([MD, (EID and ([38]) into Eq. ([22]), we obtain 

(l)a (2)a (3)a 



22]) and ([23 



Ui = u\ + u\ + u 



1 ) 



(CI) 



where 



u 



(l)a 



2E^ai 0dQ2{s) 



n 



-w 



ds 



^E^ai ^dQ2{s) 



R 



-u 



ds 



s=R 



s=R 



da'n'^n^ 



(C2) 



16 



u 



(2)a 



2E,a, . ,dQs(s) 



n 



-w 



ds 



2E^a,.^dQs{s) 



R 



-u 



ds 



--R 



s=R 



da'n"n^ 



(C3) 



u 



(3)a 



■—^u u^v? I dan ^^^^^^Q,(s) 
2Eiai 



R 



ds^ ds^ ds^ 
6 dQi 6 d'^Qi 3d^Qi 



u \u\ 



with El = 4i?/(15?7o). In these derivations, we have used the following relations; 
— / dan°'n'^n'n^ = - {5ap5^s + da-^dps + das^p-^i) , 



R 



da-rfn^ = 6as- 



(C4) 

(C5) 
(C6) 



In order to calculate U2 in Eq. ( l23l) , we need the gradient of the concentration c. 

„ 1 a 7X <9Q2 , „ ^52g2 



V^c(2) 



9sT ds^ 



-u 



ds'y ds°' dsP 

dQ3 , „ ,52g3 



-^(V-^°^")^ + ^"n^ 



ds"^ 



(C7) 



(C8) 



V7c(3) = u"u^u'J-J-^^Q,{s) 



3 



asT as" dsl^ ds^ 

n' m + 3 m u n n^ 



+ 3n^n"n^n"n^ - 5n'^n'^n''n°n^n''n^ 



5Q4 ^a2g4 



ds 



ds^ 



+ — |_|n| u n rC + v/u u'^n n^ — 2u u'^u n n^rCn \ 



R 



ds^ 



ds^ 



Substituting these into Eq. (!23|) . we obtain 



u 



(i)„ _ 2X0-1 a^Q 



R 



u 



ds 



s=R 



u 



(2)« 



2Xai.^dQ^ 



u 



R 

(3)a _ 2X0-1 

R 



-u 



ds 
I QdQi 16 d'^Qi 3 d^Q^ 



I 1 2 a 

ki u , 



i?2 5 5s ' R5 ds'^ 5 9s3 
where X has been defined by Eq. ( iBTOJ) . Comparing Eqs. O-dCl) with Eqs. (ICTOD 
we note that U2 = — (3/4)ui. 



(C9) 

(CIO) 
(Cll) 
(C12) 
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Appendix D: Correction from the convective term 



In this Appendix, we calculate the coefficient r by taking account of the correction from 
the convective term in Eq. (^. Up to the first order of V ■ {vc), Eq. flU]) has an additive 
correction as 

.dQ2 



TU 



-M- 
MR 



ds 



u 



s=R 



da'tf 



L^q 



Gqe-"^<^^+^ / d^re''^'^' {v{r')-iqc^^\?) 



(Dl) 



where we have used the relation U2 = —{3/4)ui. The vector v{f) in the second term is 
the velocity field around (and inside) the droplet moving at a constant velocity u along the 
2;-axis and is given by2^ 



v{fG + r) 



u 



u 



5 3r£\ „ I 3z „ 

2 if 2 / ^^ ^ 2R^ ^ 



-^ le + 



3R^2 

2r5 



{for r < R) 
{for r > R) . 



:d2) 



Analytical evaluation of the integrals in Eq. (IDip seems impossible in a general condition. 
Here we consider the limit R = (3R -^ 0. In this case, we may approximate G{f) as 
G{'r) = l/{4'iTDr) and r is calculated as 



r 



MR Att 
n 27^2 



AR^ 



3 /I 1 
2 U ~ 7 



1 



MR 4vr_^5^_93 



(D3) 



n 271)2 140' 

If the second term in Eq. (IDip is ignored, the factor 93/140 is replaced by 6/5. 

We can also calculate the coefficient m by taking account of the correction from the 
convective term in Eq. (jH]). 



mu 



-M- 



ds 



+ 



MR 



da'n^ 



u 

=R 



d^re^'^r' ^7(f'/).^^c(°)(P) 



(D4) 



The lowest order contribution from the ffist term is given by 



M^Q^{s) 



--R 



AMR^ 1 



( terms finite for /2 — )■ 



(D5) 



24^3 p^ V 

The second term due to the convection of the composition c has no term which is infinite for 
/3 — )■ 0. Thus, the contribution to the coefficient m from the convection of component c is 
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found to be higher order of PR. We expect the same situation for g but have not confirmed 
it since the expression is very comphcated. Finally, we make a remark that the smallness of 
e in Eq. ( I24p is independent of the smallness of R. 
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